Energy cost associated with vortex crossing in superconductors 



M.A. Moore 

Department of Theoretical Physics, University of Manchester, 
Manchester, Ml 3 9PL, United Kingdom 

N.K. Wilkin 

Department of Physics, University of Sheffield, 
Sheffield, S3 7RH, United Kingdom. 
(29 th April 1994) 

Abstract 



Starting from the Ginzburg-Landau free energy of a type II superconduc- 
tor in a magnetic field we estimate the energy associated with two vortices 
crossing. The calculations are performed by assuming that we are in a part of 
the phase diagram where the lowest Landau level approximation is valid. We 
consider only two vortices but with two markedly different sets of boundary 
conditions: on a sphere and on a plane with quasi-periodic boundary condi- 
tions. We find that the answers are very similar suggesting that the energy 
is localised to the crossing point. The crossing energy is found to be field 
and temperature dependent - with a value at the experimentally measured 
melting line of U x — 7.5kT m ~ 1.16/c^, where cl is the Lindemann melting 
criterion parameter. The crossing energy is then used with an extension of the 
Marchetti, Nelson and Cates hydrodynamic theory to suggest an explanation 
of the recent transport experiments of Safar et al. . 
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I. INTRODUCTION 



Thermal fluctuations are widely believed to play an important role in the physics of high 
temperature superconductors. In this paper we consider a possible manifestation of these 
fluctuations — that the vortex lines will be able to cross through each other. We begin by 
estimating the energy associated with such a process starting from the Ginzburg-Landau 
free energy functional. We then analyse some of the transport measurement data of Safar 
et al. 0, which has connections with the flux line crossing energy, using a modification of the 
Marchetti, Nelson and Cates hydrodynamic theory for vortex motional! 

The crossing energy calculation has previously been attempted within London theory, 
an approximation which is valid when a typical vortex separation is much greater than the 
vortex core size, and by considering just two vortices in a fixed background@~@. However, 
London theory has an unsatisfactory feature in that it requires a cut-off to be introduced at 
the length scale of the core radius. This is a difficulty as for crossing the cores of the vortices 
need to overlap. We shall instead provide an estimate of the crossing energy of two vortices 
using Ginzburg-Landau (GL) theory, and the LLL (lowest Landau level) approximation. 
Due to numerical difficulties our system is also restricted to two vortices, but with two very 
different choices of backgrounds: we have them moving on a sphere, and also on a plane with 
quasi-periodic boundary conditions. In the planar case we consider a unit cell commensurate 
with a triangular lattice ground state. We then calculate the energy difference with respect 
to the perfect vortex lattice. The results of these calculations are similar, suggesting that 
the extra energy of two crossing vortices is localized to the crossing region. 



II. REVIEW OF LANDAU-GINZBURG THEORY RESULTS 

To fix notations, we briefly describe anisotropic Landau-Ginzburg theory for a supercon- 
ductor, where xp, the wavefunction is our spatially dependent order parameter. For a fuller 
explanation and justification seeH!. 

We start with the free energy functional, 

- ' d'r (a(T) W + ?W + £ \(-m Stl -2eAMf) + * . (1) 



k B T c J \ 2 ^ 2m M 

Here a(T) is the temperature-dependent variable, (3 is the coupling constant, and is the 
effective mass. In the cases we consider the masses in the a6-plane are taken as equal and 
are denoted by m^, and the mass in the c-direction is written as m c . The temperature 
dependence of ct{T) is taken to be linear, a{T) = (T — T c )a' . We also assume the lowest 
Landau level (LLL) approximation, which is valid near the H C2 line. If we allow the variation 
in B to be determined by the equation B = HqH + {fioeh)/(m ab )(\tlj\ 2 ) (valid near to the 
H C2 boundary) then we can write our temperature variable as a B = a + efi Hh/m ab . This 
is zero along the H C2 line. A neater variable to work with is a^, which is dimensionless. It 
is related to an by, 

aH= { 4^ ^ (2) 
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The temperature dependence of the «t variable is such that high temperature is repre- 
sented by ax — > oo, low temperature by a? — » — oo and ay = corresponds to being on the 
H C2 line. The above treatment is valid where the LLL approximation can be trusted which 
probably requires at the very least that H > H C2 / f, where / = 3 according to Tesanovic et 
alM The LLL approximation allows us to write the order parameter as a set of orthonormal 
basis functions, the form of these functions being determined by the boundary conditions 
imposed to account for the finite size of the system. 



III. DETERMINATION OF THE CROSSING ENERGY 

In this section we outline the method used for calculating the crossing of two vortices 
within the GL framework, and using the LLL approximation. Furthermore our entire system 
only contains the two crossing vortices. For such a small system there is clearly a risk that the 
crossing energy will be heavily dependent on the boundary conditions used. Hence we have 
carried out the calculation with two different kinds of boundary conditions, firstly with the 
vortices confined to the surface of the sphere and then on a plane filled with unit cells with 
quasi-periodic boundary conditions. The results of our two calculations yield remarkably 
similar results — indicating that the crossing energy is relatively insensitive to the presence 
of the other vortices in the system. 



A. On a sphere 

We first consider the two vortices moving on a surface of a sphere, a geometry that 
has been shown to reduce finite size effects in numerical studies of superconductivitySIll 
and the quantum Hall effect^. The disadvantage of this approach is that the sphere is not 
commensurate with the triangular (Abrikosov) lattice ground-state. However, as we will 
show later, the energy of our ground-state is not very different from that of the triangular 
lattice. 

On a sphere we use the formalism of O'Neill and Moordll. The starting line is the free 
energy functional given in Eq. ([!]). We then place a monopole at the center of the sphere, 
which produces a radial magnetic field, satisfying BAnR 2 = A$o, where N is the number of 
vortices, which is two in our case. A choice of vector potential compatible with this condition 
is: A r = Aq = 0, = BRtem(9/2), in the usual spherical polar coordinates. The order 
parameter ip(0, <f>) can then be expanded as eigenstates of the operator (— ihV — 2eA) 2 /2m a &. 
Within the LLL the orthonormal set can be written as: 

TV 

^(0,M)= EvmCOVUM) (3) 

m=0 

where, 

i/> m {0, 0) = k m sin m (#/2) cos N ~ m (9/2)e m ^ , (4) 

with m = 0,1. ..A; k m = [(N + l)!/47ri? 2 m!(A — m)!] 1 / 2 and h is the distance along the 
vortex line and runs from — oo to oo. The zeros of ip correspond to the positions of the 
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vortices at height h. The requirement for 'crossing' is that the two zeros lie on top of each 
other at some value of h. We have calculated the energy associated with the vortices being 
in their equilibrium configuration at both ends and 'crossing' at h = 0. This has been done 
by writing down the Euler-Lagrange equations associated with two vortices and then solving 
with boundary conditions at h = and h = oo compatible with the 'crossing' requirement. 



(Technical details of the calculation can be found in Appendix [A~l 



To follow the above procedure we need to write the free energy in terms of the coefficients 
v m . The general expression becomes@, 



F[{Vm}] 



dh 



f3B 



E( 

. m=0 

N 



h 2 , dv. 



2m r dh 



m |2 i I |2\ 



J2 W(m + p,m, n)v m v p v* n v*8, 



m+p,n+r 



m,n,p,r=0 



(5) 



where 



W(m + p, m, n) 
and f(x,y,z) 



2(N + l) 2 f(m + p, m, n) f(2N — m — p, N — m, N — n) 
N(2N + 1) f(2N, N, N) 

x\/[(y\z\(x - y)\(x - z)\) 1/2 2 x+2 } . 



(6) 



From Eq. (|5|) and Eq. (^|) we find that the free energy/per unit length in the z-direction for 
two vortices can be written as: 



dh 



dv 



* |2 



i=0 



2m r dh 



+a H | Vi | 2 ) + /35/$o(0.45M 4 + °- 3 kil 



+ 0.45H 4 + 0.9|w |>i| 2 + 0.3|w |>2| 2 + 0.9M 2 H 
+0.3(v o v 2 v*vl + vlvlvyvx)) 



(7) 



The first task in evaluating the crossing energy AF, is to find the baseline, the free 
energy minimum of the ground state of the system. It is easily found by minimisation of 
Eq. (0) to be — (ot 2 H <&o) / {1.2(3 B) per unit length, as compared with the Abrikosov ground- 
state of — (a 2 1 §o)/(f3 a f3B). Hence we have an effective f3 a = 1.2, which is a reasonable 
approximation to the Abrikosov triangular lattice value of 1.16. This solution corresponds 
=l=c, where c = 



to Vo = ±c, v 2 



Offl^o/l-2/3-B. To find the crossing energy we next 



write down the Euler-Lagrange equations for the coefficients Vi(h). Solutions exist for which 
v i = and in which Vq(— oo) = — c, v (oo) = c, f 2 (— oo) = c, 1*2(00) = c. The crossing point 
then occurs at h — when Vq — 0, which happens in our notation to be when both vortices 
are at the south pole of the sphere. (See Appendix |A 1| for details.) 

The idea was to start from the crossing point and integrate forward numerically along h 
to the ground state configurations. Using energy conservation and symmetry arguments at 
h = left us with one remaining unknown parameter there. We then varied this parameter 
until we found a value that allowed us to integrate forward and obtain our equilibrium 
configuration. Classically this can be thought of as finding the trajectory of a particle that 
is stationary at the maximum in the potentials at (-c,c) to (c,c) , which moves as dictated 



4 



by the Euler-Lagrange equations. (See Fig. (P for this potential, and the actual trajectory- 
taken.) 

Having found a suitable trajectory, we could numerically integrate along the path to find 
the difference in energy from the equilibrium state. This left us with a prediction for the 
crossing energy of: 

1.3 I a H | 3/2 ft$ . ,, /9 , m 

AF = ' H = 0.33 a T 3/2 kT. 8 

Our procedure could be generalized to handle several vortices. For two vortices, use of 
symmetry and conservation arguments resulted in only one adjustable parameter, but with 
three vortices (say) we would have at least a three dimensional space to search. The problem 
rapidly becomes completely intractable as the number of vortices increases and hence we 
did not attempt to go beyond two. 



B. On a plane 

Given that we could not easily solve for more than two vortices on a sphere, we instead 
looked at the same problem under very different boundary conditions. Fortunately we obtain 
a very similar answer which suggests that the crossing energy is determined mainly by the two 
vortices which are being crossed. This second set-up was on the plane, using the boundary 
conditions suggested by Kato and NagaosaQ. A rectangular unit cell is chosen with sides 
commensurate with the formation of a triangular lattice of vortices, in effect forcing the 
ground-state to be a triangular lattice. The cell also has quasi-periodic boundary conditions 
such that the same motion of the vortices is carried out in all of the unit cells simultaneously, 
which at first sight looks rather unrealistic! 

Within the LLL, and working in the gauge (0, Bx) the order parameter can be written 

as, 

AT-l 

^(x, y,h) = J2 c n (h)(j) n (x, y) , (9) 

n=0 



where the <p n are defined by (see Ref. 0): 

\ *. 1 / \ 2 ~ 

(10) 



/ 1 

<Pn{x,y) = J- — j=- ex P 

V L V V 7 ™ m=-oo 



'27T/ 2 , , \ * l( 27T/ 2 



and I 2 = L x L y /(47r). L x and L y are the sides of the cell. The ground state is a triangular 
lattice if y/3L y = L x . The periodicity of n is such that 

L x y 

<p{x,y + L y ) = <p(x,y) and <j>{x + L x , y) = <j>(x, y) exp(-z-^-). (11) 

If we substitute the order parameter into Eq. ([[]), for general iV (number of vortices in the 
system) we obtain an expression for F[{c n }]. We write this using some auxiliary notation, 
explained below; 
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roc N ~t 

F[{c n }] = dhJ2 

J -°° n=0 



2m r 



den 
dh 



+ ot H \c n \ + 



^\hK 1 Tc c c * c * e x P r-^(p 2 + Q 2 

4$oV L y (2M + l) 2 f CniCn2V ^ eXP { 2N [r +W 



(12) 



The sum indicated by the symbol Y^,i runs over the eight variables ni, . . . n 4 and mi, . . . 1714, 
subject to the requirements that 



n 1 +n 2 + N(mi + m 2 ) = n 3 + n 4 + N(m 3 + m 4 ) 



(13) 



< n < N - 1 
The symbols P and Q are defined as 



and - M < m < M 



P = rii — n 2 + N{m\ — m 2 ) 
Q = n 3 - n 4 + N(m s - m 4 ) . 



(14) 



(15) 
(16) 



If we now consider only two vortices and evaluate the sums, we are left with the equation 
for the free energy in the form: 

i=0 



+q\c \ 2 \ Cl \ 2 + s(c 2 cl +c* ci)) 



(17) 



where p = 0.465337, q = 0.245010 and s = 0.008063. The ground state energy per unit length 
is F = — ajj^o/ j3 a /3B with (3 a = 1.15956. This is the triangular lattice result, accurate to 
five decimal places, the error occurring in the evaluation of the sums in <f> n . 

We now outline the calculation required to evaluate the crossing energy on the plane, 
leaving the details for Appendix [A 2j . Symmetry suggests that a suitable place for the zeros 
of the vortices to be superposed is at the center of the unit cell. To achieve this we imposed 
the condition, 



Q>(O)0o ('y,^)+ci(O)0 1 (y,^)=0. 



(18) 



The equilibrium configuration can be satisfied with co( 
■00) 



-00 



C, Ci( 



-00 



-c, Col 



-00 



Ci{— 00) = c where we know c = y\ajj\^r,/ (3 a (3B by energy arguments. By imposing 'energy 
conservation', and recognising that we also have 'angular momentum conservation', in the 
dynamical analogue of a particle moving in a potential it is possible to reduce our set of 
Euler-Lagrange equations again to solving for just one unknown parameter. We can then 
proceed in a similar manner to the case of the sphere, although in this case we have numerical 
problems in reaching the equilibrium condition and the energy we evaluate has a residual 
kinetic energy which is 5% of the maximal kinetic and 2% of the total energy, and hence an 
error of ~5% on the crossing energy. The final value we calculate is very similar to that for 
the sphere, Eq. (H), 
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(19) 



It should be remembered when comparing the two calculations of the crossing energy 
that we already had a 4% difference in the ground-state energies for the different boundary 
conditions. 



IV. CROSSING ANGLE 



Previous calculations of the crossing energy within the London regime by Obukhov and 
Rubinsteinlll and Nelsoni have considered a mechanism for crossing in which at the crossing 
point the the configuration of the vortices is chosen for simplicity to be that shown in Fig. E3. 



The crossing angle 6 C depends on the anisotropy of the sample, 6 C ~ arctan(ym c /m a fc), and 
is 45° for an isotropic superconductor. 

In our calculation it was not necessary to impose the configuration of the vortices at the 
crossing point, but in the planar case it can be extracted from the solution to the Euler- 
Lagrange equations. (For the case of the sphere it is difficult to decide how the angles map 
to those of the flux lines physically crossing.) It should be noted that the lowest energy 
configuration for the vortices to cross may differ between the London and LLL regimes as 
in the latter the variation of the -B-field is not included. 

On the plane, we know the positions of the zeros at all heights is given by, 

c (h)MX, Y) + Cl (h)M x > Y) = 0, (20) 

where co(h) and c\(h) are known from solving the Euler-Lagrange equations. In order 
to investigate the behaviour of the flux lines near the crossing point we consider deviations 
(x,y) from the position at h — 0, X — L X /2,Y = L y /2, by considering the differential of 
Eq. (pO"D with respect to h; 



{d h c )<f)o + (4ci)0i +{d h x) {c d x (j) + cA^i) +{d h y) {c d y <f) + c^yfa) = 0. (21) 

We find that at the crossing point a = 0.565i and that (3 = 5.26(x—iy) and 7 = — 5.26(ix+y). 
Hence solving for the two equations for the real and imaginary parts we find that in the 
vicinity of the crossing point x = y = y/O.llh. 

This implies that the vortices save energy on close approach by allowing their core en- 
ergies to cancel (see Fig. ||). The vortices also approach not along the body diagonal of 
the unit cell but at 45° such that there is a twist along the length of the vortices as they 
approach. 



V. COMPARISON WITH LONDON THEORY RESULTS 

In order to allow comparison between our results for the crossing energy and the estimates 
for crossing within the London regime, it is simplest to re-express our crossing energy in 
terms of the Lindemann melting criteria. This criteria suggests that when the thermally 
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induced spatial fluctuations in the positions of the flux lines (u 2 h ) become of the order of 
c^Oq the lattice will melt. The Lindemann number can be obtained indirectly from neutron 
scattering and estimates range from 0.1-0.4. 

By writing the Lindemann criterion in the form suggested by MooreS, 

M ^^^ = clal a » = l (22) 

and using the formulae for 

\a H \h 2 0.2A\a H \ 2 

cm = — — > (23) 



fS m c (3(3a 00 13113 

which are determined using elasticity theory,0 we find that c\ = 3.13/|ar| 3 / 2 for |ar| at the 
melting line. Hence, using the result for U x ~ 0.35|a r | 3/2 fcT, we find that U x /kT ~ 1.1/4 
(leading us to estimate cl ~ 0.4). Nelsoni has estimated in the London regime that 



i^ 2( V2-l).& ,n M (*) (24) 
k B T y m c \4:Tc\ ab J c£ 

which is of the same order of magnitude as our LLL estimate. A more recent estimate, 
also in the London regime, by Carraro and FisherB leads to a result of U x /kT * 0.24/4, 
which is rather lower than our result. There is of course no reason to expect that the results 
for the crossing energy in the LLL and the London regime to be identical. 



VI. DISCUSSION 

Having estimated the crossing energy we now consider its relevance within the flux liquid 
phase of the high temperature superconductors. In order to understand the magnitude of 
the crossing energy in terms of physical quantities we first estimate its value in terms of 
fcTjvf at the melting line found experimentally by Worthington et a/.@ and Safar et a/.0. 
(Providing we accept that the melting line can be equated with the irreversibility line in 
magnetisation measurements.) The latest data from Safar et al. has been taken in fields up 
to 16 T. Their results show a first order transition for fields below ~ 10T and a continuous 
transition (which is no longer a melting transition) for larger fields. In order to consider the 
magnitude of the crossing energy in the vicinity of the transition we need an estimate of a^. 
From previous workS we have a crude estimate of «t ~ —8 based on the melting line shown 
in the Worthington et a/.0 data. The latter data is taken in the region now associated with 
the first order transition but the value of olt should not change appreciably with increasing 
field. An alternative estimated from the analysis of theoretical specific heat data suggests 
ay ~ — 7, although their argument is somewhat suspectfi If we believe in this magnitude 
of Q!t then we have U x /kT ~ 7.5 (for \oct\ = 8) in the region of the phase diagram being 
investigated. 

Previously Marchetti, Nelson and have shown using a hydrodynamic theory 

that the effects of a twin boundary or similar pinning surface whose normal is perpendicular 
to the field act over a characteristic length S a b, when the magnetic field is in the c-direction. 
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This characteristic length is effectively the distance between line crossing events in the ab- 
plane. CatesS found that 5 a b — a exp[U x /2kT], where a is the spacing between vortex 
lines. By extending Cates argument we now deduce the length scale S c , relevant for velocity 
gradients varying along the field direction, which we believe to be the important length 
scale in the transport measurements of Safar et al. 0. The length 5 a b is the typical transverse 
distance between crossing events, and can be written as 



lp 



where this is just a 'random walk' of 8 c /l p steps. l p is a persistence length, as described 
by Cates, and is related to l e , the distance along a flux line which one has to travel before 
encountering another flux line by l e l p = a§. The resulting distance between crossing events 
in the c-direction is now given by: 

6 c ^l e e u «/ kT . (26) 

It is clear that when 5 C becomes equal to the thickness of the sample (or larger) then the 
vortex motion at the top face of the sample will be strongly correlated with that at the 
bottom face of the sample. Experimentally the degree of correlation between the motion of 
vortices at the top and bottom surfaces of the sample has been measured by Safar et al. on 
YBCO crystals and has been analysed by Huse and Majumdaril. The experimental set-up 
involves applying a magnetic field along the c-axis and then injecting a transport current 
in the top a-b plane along the a-axis. The voltages in the top and bottom faces are then 
measured and the maximum field temperature combinations above which V top = V^ot are 
recorded. If the vortices are not readily able to cross through each other then the voltages 
on the two faces due to the Lorentz force induced on the top surface should be the same. 
However, once 5 C = L, the size of the sample, line crossing will allow the vortices to remain 
pinned on the bottom face whilst they continue to move due to the Lorentz force on the 
top face. The mechanism for this procedure can be seen in Fig. Initially, the flux line 
labeled 'A-B' is in front of the line 'C-D' and the top of line 'A' is subject to the Lorentz 
force exerted by the transport current. However, line 'C-D' is pinned along its entire length 
and is unmoved by the Lorentz force. If crossing is energetically favourable then when end 
'A' encounters end 'C the two lines will cut and re-combine as 'A-D' and 'C-B' resulting in 
a net movement of vortices, and hence a voltage in the top surface. There is no net flux line 
movement and hence no voltage in the bottom surface. 

Using the results from Safar et al. and our value of the crossing energy we find that 
S c = 2.5l e exp[U x /kT] is a good fit to the data for the L = 30/im sample. The value for 
l e was estimated by assuming that the persistence length l p would be of the order of the 
spacing between the Cu-0 planes in YBCO. Considering that l e is only an order of magnitude 
estimate, this mechanism is a possible explanation of the results. We also find that for a 
given magnetic field the point at which V top = V^ot decreases, with increasing thickness 
of sample. This is in qualitative agreement with experimental data at 1 T. However, the 
quantitative comparison is poor, which is not surprising as 1 T is outside the expected region 
of validity of the LLL approximation. 

Clearly, this is a first estimate of the crossing energy, and a more complete calculation 
will consider a larger number of vortices, as well as allowing for fluctuations of the vortices. 
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This calculation has also indicated that the length scales in the problem are of the order of 
the Cu-0 layer spacing indicating that an investigation of crossing within a layered structure 
is called for, rather than as here just within the continuum GL theory. 

We would like to thank Hugo Safar for the use of his experimental data and Michael 
Cates and William Barford for many useful discussions. 
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APPENDIX A: DETAILS OF THE CROSSING ENERGY CALCULATIONS 

In this appendix we give a more detailed description of the procedure used to calculate 
the crossing energy on both the sphere and the plane. The problem we had to solve was 
essentially a classical mechanics problem of finding a trajectory for a ball to roll from one 
metastable position on one local maximum to another at the same height, with a set of 
conditions to control the velocity and position at the central point of the trajectory. 

1. On a sphere 

The first simplification we used on working on the sphere was to invoke a reduced set of 
the variables {v m }. We claimed in § [III A| that this subset was sufficient to yield the correct 
solutions — the justification is quite simple although algebraically a little heavy. 

We begin with the free energy in terms of the full set of {v m }, as in Eq. (J|), and the 
order parameter if) as defined in Eq. (|3]). We then express if) in terms of spinor variables. In 
the case of two vortices, N=2 this is simply: 

if)(9, (ft) = A{wu\ — uwi){wu2 — uw 2 ) , (Al) 

where u = cos(6 l /2)e _ ^ //2 ,w = sm(9/2)e t( ^^ 2 and A = re 1 ? is a complex variable. The 
positions of the vortices are described by the zeros of if> which occur when u = u% and 
w = Wi. Hence we now know the positions of the vortices in terms of the polar angles 9i and 
4>i. By comparison of Eq. and Eq. ( |A1| ) we can find the relationship between the {v j} and 
A, {ui, Wi} such that we can write Eq. (|5|) in terms of the spinor variables. A more coordinate- 
independent variable than the individual positions of the vortices is the separation between 
them, a measure of which is the scalar product n of their positions. In terms of the spinor 
variables n can be written as n = cos(6' 1 ) cos(9 2 ) + cos(0i — (f> 2 ) sin(^x) sm(9 2 ). Inspection of 
the potential energy written in terms of the spinor variables shows it can easily be re-written 
in terms of n and r: 

V = ^-(3 + n) + -^—(39 + 30n + 3n 2 ) 
4 lot) 

The kinetic energy is somewhat more complicated. However the terms that cannot be 
written in terms of n and r can be written as sums of squares. Hence they are effectively 
conserved momenta, whose contribution will be zero when we look for the minimum energy 
solution. The remaining relevant terms in the kinetic energy are then: 

T = ^r^(d h r) 2 + - (d h r) (d h n) + —f- (d h nf. (A3) 
4 4 lo(l + n) 

If we then derive the Euler-Lagrange equations for this new system we find they map back 
to just the real parts of vq and v 2 - Hence the set of variables we have chosen to work with is 
sufficient. We should perhaps point out that there is good reason for working with the less 
intuitive variables vq and v 2 — the Euler-Lagrange equations are then far simpler. 

We are now in a position to calculate the crossing energy. We start by writing down the 
Euler-Lagrange equations, 
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(A2) 



h 2 ( 3B 

^ d l v o = v i-\a H \ (0A5v 2 + 0.3v 2 ] 



11 -d 2 h v 2 =vJ-\a H \ +^ (0A5v 2 + 0.3v 2 )) . (A4) 



2m c a \ $ c 

Then the total energy, T + V = $o|«/f| 2 /(l-2/9-B), does not vary with h. This is an initial 
value problem so we need to know the set of {^0,^2,^/1^0,^/1^2} at h=0, where the vortices 
cross. Examination of the boundary conditions showed that vq was antisymmetric and t>2 
symmetric about h=0 (see Fig. ||) and as a consequence v = and duv 2 = at h = 0. Then, 
by conservation of energy, we knew that dhV\ is a function of v%. Hence we just searched 
the one parameter space vo until we found a solution of the equations whose end-point was 
compatible with the ground-state configuration. The value of h at which the vortices became 
straight again is h ~ 6h/^2m c aH- That is the crossing was completed on a length scale l c 
of the order of the conventional phase correlation length £ c , and l c ~ 12£ c , see Fig. [| In 
practice this length scale will be of the order of the spacing of the Cu-0 planes, indicating 
the need to go beyond the continuum GL approach of this paper. 



2. On a plane 

The procedure for calculating the crossing energy on a plane is very similar to that just 
described for a sphere. Hence we will consider it only briefly. In this case no attempt was 
made to reduce the set of variables, so we had as our starting point a set of four coupled 
non-linear Euler-Lagrange equations, 

d 2 h r = r (2pr 2 +r\{q + 2s cos(20 o - 200) + (40 o ) 2 - l) (A5) 

d 2 h n = n {2 P r\ + r 2 (q + 2s cos(20 o - 200) + (d^) 2 - l) (A6) 

r <9^o = -2sr r 2 sin(20 o - 200 - 2(d h r ) {d h <f> ) ( A7) 

n^0! = 2snr 2 sin(20 o - 200 " 2(9^0(^00 , ( A §) 

where cq = r^e 1 ^ etc., and p, q, s are as defined for Eq. fllTD - Written in this form, we can see 
by inspection, that as well as conservation of energy we also have conservation of 'angular 
momentum' as dhir^dh^o + r i<9/i0O = 0- Moreover, because the boundary conditions at ±00 
require there to be no kinetic energy, the angular momentum must be equal to zero. So 
far we have a general set of equations describing the motion of the two vortices within the 
unit cell, but we need a further constraint to make the vortices cross. Having decided that 
a suitable place for the vortices to meet would be the center of the unit cell, we find that 
we need to satisfy the condition, c o (O)-0o(-^x/2, L y /2) + ci(O)-0i(L x /2, L y /2) = 0. Evaluation 
of the -0's shows that we require r\ = 0.13165r and 0o = 0i = at h = 0. In fact the 
other obvious places for the vortices to cross, (±L x /2,0) and (0, ±L y /2) yield equivalent 
constraints. Through symmetry arguments we find that c^co = dhC\ = at h = 0. 

Using all of these facts leaves us again with only one free parameter in our initial con- 
ditions. Finding a solution which reached the necessary end-point proved more tricky this 
time, and we were unable to produce a solution which bettered having 5% of the maximal 
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kinetic energy (2% of the total energy) in 'equilibrium', and correspondingly the equilib- 
rium separation of the vortices is only accurate to 2%. Hence an estimate of the error on 
the crossing energy of 5% would seem reasonable. 
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FIGURES 



FIG. 1. Contour plot of the 'potential' as described in Eqn. ( |A^ ) for the spherical boundary 
conditions problem. The solution to the Euler-Lagrange equations which corresponds to the vortices 
crossing is that linking the points 'A' and 'B', and c is as defined in the text. 

FIG. 2. Ansatz configuration for vortices crossing in the London regime used in Ref. ^, where 
7 = (m c /m ab ) 1/2 . 

FIG. 3. Configuration for vortices crossing within the LLL regime, l c is the length over which 
the crossing takes place. 

FIG. 4. Possible vortex motion in the Safar et al. experiment. Initially (a) line 'A-B' is in 
front of C-D and subject to a Lorentz force at end 'A', and line 'C-D' is pinned. Crossing enables 
end 'A' to carry on moving, (b), resulting in a net voltage difference between top and bottom of 
the sample. 

FIG. 5. The functions Vo(h) (antisymmetric)and V2(h) (symmetric) as calculated by solving the 
Euler-Lagrange equations, vq and V2 are in units of 'c' and h is in units of h/^2m c aH, the c-axis 
correlation length. 
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